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Abstract 

The cross section for the production of Z boson pairs is measured using the data 
collected by the L3 detector at LEP in 1999 in e + e~ collisions at centre-of-mass 
energies ranging from 192 GeV up to 202 GeV. Events in all the visible final states 
are selected, measuring the cross section of this process. The special case of final 
states containing b quarks is also investigated. All results are in agreement with 
the Standard Model predictions. 
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1 Introduction 



The increase of the LEP centre-of-mass energy, a/s, beyond the Z pole has extended the range 
of the accessible physics processes to include a sizable fraction of four-fermion events. An 
important part of the four-fermion final states emerges from the pair production of W or Z 
gauge bosons. 

The study of Z boson pair-production is of interest as it offers a further test of the Standard 
Model of the electroweak interactions [IL in the neutral gauge boson sector. Moreover, this 
process constitutes a background in the search of the Standard Model Higgs boson. In addition, 
Z-pair events allow the investigation of possible triple neutral gauge boson couplings, ZZZ and 
ZZ7 PP] , forbidden at tree level in the Standard Model. These events can also test new theories 
beyond the Standard Model such as Supersymmetry [Hlllj or extra space dimensions [Sj. 

At the lowest order, Z pair-production proceeds via two t-channel Feynman diagrams with 
an internal electron leg. Considering the Z decays into fermions, this process is conventionally 
denoted as NC02, from the acronym of the neutral-current production mechanism of the four- 
fermions and the number of diagrams. A wider definition is used in this letter, encompassing the 
regions of the full four-fermion phase space compatible with the pair-production of Z bosons. 
Results in the NC02 framework are also given. 

The experimental investigation of the Z pair-production is made difficult by its rather low 
cross section, compared with competing two- and four-fermion processes, that constitute large 
and sometimes irreducible backgrounds. This process was observed at threshold by the L3 
collaboration at y/s = 183 GeV [Oj and studied later with a higher statistical sample at y/s = 
189 GeV ;7J. Results from the other LEP collaborations were also reported |H1[H|- This letter 
describes the extension of the L3 analyses to centre-of-mass energies between 192 GeV and 
202 GeV. The measurement of the cross section is presented together with other results that 
include lower centre-of-mass energies. The measurement of the cross section for the particular 
case of Z pair-production and decay into at least a b quark pair is also discussed. 

2 Data and Monte Carlo Samples 

The data under investigation were collected in 1999 by the L3 detector pH] at four different 
centre-of-mass energies, 191.6 GeV, 195.5 GeV, 199.5 GeV and 201.7 GeV with corresponding 
integrated luminosities of 29.7 pb _1 , 83.7 pb _1 , 82.8 pb -1 and 37.0 pb _1 . These energies are 
denoted as 192, 196, 200 and 202 hereafter. 

The EXCALIBUR [H] Monte Carlo is used to generate events belonging to both the sig- 
nal and the background neutral-current four-fermion processes. Background from fermion-pair 
production is described making use of PYTHIA [T2j and KK2f (THj (e + e _ — > qq(7)), KO- 
RALZ [Hj and KK2f (e+e~ -> and e+e" -> T + T~(i)) and BHWIDE P3j (e+e" -> 

e + e~(7)). Background from charged-current four-fermion processes is generated with EXCAL- 
IBUR for ez/ e qq' and £ + Pi£~&e with £ = e,fi,T and KORALW jT^] for W pair-production and 
decay in the final states not covered by the simulations listed above. Contributions from 
multi-peripheral processes are modelled by PHOJET [T^j (e + e~ — > e + e~qq) and DIAG36 [IE] 
(e + e~ — ► e + e~£ + £~), in the quark and lepton low invariant mass region not included in the 
samples generated with EXCALIBUR. 

The L3 detector response is simulated using the GEANT program which takes into 
account the effects of energy loss, multiple scattering and showering in the detector. Time 
dependent detector inefficiencies, as monitored during the data taking period, are reproduced 
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in these simulations. 

The Z pair-production signal is defined as the subset of the full four-fermion phase space 
satisfying the following requirements [HIE!- First, the invariant mass of both fermion pairs 
must be between 70 GeV and 105 GeV. In the case in which fermion pairs can originate 
from a charged- current process (uddu, esse and u^£ + i7i£~, with I = e,/i,r), the masses of the 
fermion pairs which can also emerge from W decays are required to be either below 75 GeV or 
above 85 GeV. Finally, events with electrons in the final state are rejected if for any electron 
| cos# e | > 0.95, where 9 e is the electron polar angle. 

The expected cross sections for the different final states are computed imposing the require- 
ments described above on a sample of events generated with EXCALIBUR, and found to be 
0.79 pb, 0.92 pb, 1.00 pb and 1.03 pb at the four centre-of-mass energies, in increasing order. In 
this calculation a s = 0.119 j2H| is included for the QCD vertex corrections. An uncertainty 
of ±2% is assigned to these predictions, reflecting the differences between them and those ob- 
tained with the GRC4F [21] Monte Carlo generator as well as the expected accuracy of the 
treatment of initial state radiation. 

The relative populations of the different channels, as obtained from their corresponding 
cross sections, differ slightly from those of the NC02 framework, derived from the Z branching 
ratios [20]. 

The cross section for final states with at least one b quark pair is significantly smaller than 
the total cross section and all centre-of-mass energies are hence combined. The corresponding 
predicted cross section is 0.27 pb with an uncertainty of ±2%. 

3 Event Selection 

All the visible final states of the Z-pair decay are investigated with criteria similar to those 
used at y/s = 189 GeV |7], and modified to follow the evolution of the signal topology, which 
manifests a larger boost of the Z bosons with the higher y/s. All selections are based on the 
identification of two fermion pairs, each with a mass close to the Z boson mass. 

Electrons are recognised from energy depositions in the electromagnetic calorimeter whose 
shower shape is compatible with that initiated by an electron or a photon. According to 
the selection channel, a track as reconstructed in the central tracker may be required to be 
associated to this cluster. 

Muons are reconstructed either from tracks in the muon spectrometer pointing to the in- 
teraction vertex and in time with the event, or via only energy depositions in the calorimeters 
consistent with a minimum ionising particle (MIP) which have a matching track in the central 
tracker. 

Tau leptons are identified from their decay into electrons or muons or as low-multiplicity 
jets with one or three associated tracks. A total unit charge is required. 

Quarks manifest themselves with a high multiplicity of calorimetric clusters and charged 
tracks. These are grouped into the required number of jets by means of the DURHAM algo- 
rithm [221. 

The hermeticity of the detector allows to reconstruct the four-momentum of Z bosons de- 
caying into neutrinos by means of the event missing energy and the momentum imbalance. 

The selection criteria specific for each final state are discussed below, first for the channels 
containing hadrons, then for the purely leptonic ones. 
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3.1 The qq£ + £~ Channel 

For each of the final states qqe + e~, qq/i + /U~ and qqr + r~, a dedicated selection is performed. 
A pair of leptons should be present in high multiplicity events with visible energy and effective 
centre-of-mass energy respectively in excess of 0.5-y/i and 0.6y/s. The effective centre-of-mass 
energy is the energy at which the e + e~ interaction takes place after the possible emission of 
initial state radiation photons. It is reconstructed taking into account both photons observed 
in the detector and those collinear with the beam axis [23] • For the qqe + e~ selection, at least 
one electron is required to have a matched track. No more than one MIP is allowed in the 
qqfi + [i~ selection. 

The qqr + r~ selection relies on both a particle-based and a jet-based approach. The former 
is aimed to identify a pair of taus in the events while in the latter the event is constrained 
into four-jets. Two of the jets must have less than four tracks and are considered as the tau 
candidates. At least one of them has also to coincide with an identified tau. The radiative qq(7) 
background is further suppressed by rejecting events containing a photon of energy larger than 
40 GeV. 

The topology of the pair-production of Z bosons is enforced by requiring the lepton pair 
and the jet pair to have an opening angle of at least 110° for the electron and muon channels 
and 120° for the tau channel. Moreover, the invariant mass of the jet-jet and the lepton-lepton 
systems after performing a kinematic fit, which imposes energy and momentum conservation, 
must be between 70 GeV and 120 GeV, as depicted in Figure for the lepton case. 

The contribution from semileptonic decays of W pairs is reduced by requiring the transverse 
missing momentum to be lower than 0.2^/i and the visible energy in the electron and muon 
channels to be at least 0.85\/i, while it has to be between 0.6-y/i and 0.9-y/i for the tau selections. 

To reject the residual background from W pair-production and hadronic events with gluon 
radiation, the events are subject to the DURHAM algorithm requiring y 34 to be greater than 
0.001 for the electron and muon channels and 0.0025 for the tau channel. 2/34 is the DURHAM 
resolution parameter for which events change from a three-jet into a four-jet topology. 

The kinematic fit is repeated on events that pass at least one of the three selections described 
above, with the additional constraint of equal invariant masses for the jet-jet and lepton-lepton 
systems. The distribution of this invariant mass, M 5 c, is shown in Figure QJd. Table 1 sum- 
marises the efficiencies achieved by the different selections. Tables 2 and 3 present respectively 
the total yield of the selection and its breakdown into the different centre-of-mass energies. 



Final State 


Selection 


e + e qq 


/!+ /j qq 


t + t qq 


Total 


e + e~qq 
/i+/i~qq 
r + T~qq 


77.2% 
0.6% 


53.7% 
0.2% 


2.6% 
6.2% 
28.2% 


79.5% 
59.2% 
28.7% 



Table 1: Efficiency of the qq£ + £ selections and of their combination. 



3.2 The qqvv Channel 

The selection of the qqvv channel proceeds from high multiplicity events with an invariant 
mass in excess of 50 GeV. These criteria deplete the total data sample of purely leptonic two- 
fermion final states and products of two-photon interactions. Hadronic events from qq(7) and 
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Table 2: Data, signal and background Monte Carlo events selected by each analysis and their 
efficiency. The qqvv and qqq'q' entries are reported for selection requirements of 0.5 and 0.2 
on the neural network outputs, respectively. The £ + £~i/u entries refer only to electrons and 
muons. The uncertainties shown are from Monte Carlo statistics. 



W pair-production are then reduced by requiring the invariant mass to be less than 130 GeV 
and the mass recoiling against the hadronic system to exceed 50 GeV. Semileptonic decays of 
W pairs are suppressed by rejecting events with electrons or muons with energies above 20 GeV. 
The missing energy signature of a Z boson decaying into two neutrinos is further exploited by 
requiring the transverse momentum to be above 5 GeV, the energy deposition in the forward 
calorimeters to be below 10 GeV and the missing momentum vector to point at least 16° away 
from the beam axis. Moreover, the energy in a 25° azimuthal sector around the missing energy 
direction, E25, must not exceed 30 GeV. 

The selection requirements described above select 407 events in the full data sample. The 
Monte Carlo expectations are 45 events for the signal and 339 for the background, mainly 
accounted for by charged-current four-fermion processes. 

An artificial neural network is then designed to further discriminate Z pair events from 
background. It is based on event shape variables that differentiate the two-jet from the three- 
jet topology, on the sum of invariant and missing masses, on the masses of the two jets into 
which the event can be forced, the total missing momentum and E25. A constrained fit is 
applied to the hadronic system in the hypothesis that the missing energy and momenta are 
due to a Z boson. The resulting mass Mfu, presented in Figure Ct, is also used in the neural 
network. The output NNout of the neural network is presented in Figure CJl. The efficiency and 
the results of this selection are summarised in Table 2 and detailed in Table 3 for the different 
centre-of-mass energies for a benchmark cut of 0.5 on NNo u t- 

3.3 The qqq'q' Channel 

The qqq'q' channel is investigated by first selecting high-multiplicity events with a visible energy 
between 0.6-y/i and 1.4-^/i, parallel and perpendicular imbalances below 0.3\/s and no identified 
electron, muon or photon with energy above 65 GeV. The events are forced to four jets and 
then subjected to a constrained fit which rescales the jets to balance momentum while imposing 
energy conservation, greatly reducing the dependence on the calorimeter energy scale. 

Such hadronic events are copiously produced in QCD processes and W-pair production. 
Two artificial neural networks are sequentially constructed to isolate the Z pair signal and 
reject these two backgrounds. The first neural network j21j helps in selecting the signal from 
the QCD background. After a cut at 0.65 on the output of this neural network, displayed in 
Figure EK, the content in W and Z pair production events is enhanced. 

A second neural network is then built to distinguish Z pairs from W pairs. It relies on 
the reconstructed di-jet mass, the maximum and minimum jet energy, the average number of 
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Table 3: Number of data (N D ), signal (Ns) and background (N B ) Monte Carlo events selected 
at the different centre-of-mass energies in the separate final states and their sum. The qqz/z/ and 
qqq'q' entries are reported for selection requirements on the neural network outputs of 0.5 and 
0.2, respectively. Monte Carlo statistical uncertainties are given on the signal and background 
expectations. 



charged tracks per jet and the di-jet mass difference. 

Almost 40% of the events generated in the qqq'q' channel and satisfying the signal definition 
contain at least a b quark pair, while the b-content in W pair events is negligible. A b-tag 
discriminant is then added to the network to further discriminate Z pair from W pair 
events. 

Figure Eb displays the output of this network after the W pair enriched region below 0.2 is 
discarded. Events compatible with Z pair-production preferentially populate the region between 
0.6 and 0.8 if their content in b quarks is low and lie above otherwise. The performances of 
this analysis are summarised in Tables 2 and 3. 

These results are confirmed by a simpler cut-based analysis that relies on the signature of 
the different boost of Z and W pairs as retained in the two di-jet opening angles, the di-jet 
mass difference and the di-jet mean mass. Another study mainly aimed at the rejection of the 
QCD background and the simultaneous selection of four-jet W and Z pair events also yields 
compatible results. Both these analyses are affected by a lower purity which follows from the 
absence of a b-tag. 

3.4 The l + t-vv Channel 

A pair of identified electrons or muons constitutes the core of the £ + £~vv selection. Tracks are 
not required for the electron identification and MIPs are not considered as muon candidates. 
The lepton pair must be consistent with a Z boson, with an invariant mass, Mm, between 
80 GeV and 100 GeV. The recoil mass, M rec , is required to lie in the same interval to enforce 
the signature of the second Z decaying into two neutrinos. 
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Fermion-pair events are rejected by requiring the lepton pair to be acoplanar and to have a 
visible energy compatible with the signal hypothesis. Moreover, the missing momentum vector 
must point away from the beam line. 

The background from other resonant and non-resonant four-fermion processes is reduced by 
performing a kinematic fit which imposes the Z mass to the visible pair of leptons. The recoil 
mass M/g* is recalculated and required to be compatible with the Z mass. 

The distribution of the sum of Mm and M rec is expected to peak around twice the Z mass 
for signal events, and is presented in Figure Eh- The efficiency of the selection is reported in 
Table 2, which also lists the total number of selected and expected events, detailed in Table 3 for 
the different centre-of-mass energies. No contribution from the t + t~i/v signal is expected after 
this selection. The dominant background arises from charged-current four-fermion processes. 

3.5 The £+£-£'+£'- Channel 

To achieve a high efficiency, the selection of the £ + £~£' + £'~ channel starts from four or more 
loosely identified leptons in low-multiplicity events and concentrates on the kinematic properties 
of just a pair of them. Electrons with or without a matched track, muons and taus are accepted 
in the first stage, provided their energy exceeds 3 GeV. If more than four leptons are present, 
the four most compatible with energy and momentum conservation are selected. 

The event must contain at least an electron or a muon pair. To form such pairs at least 
one electron should have a matched track and no MIPs are considered as muons. In the case 
of multiple choices, the pair with the invariant mass M« closest to the Z mass is studied. Both 
Mju and the recoil mass M rec to this selected lepton pair are required to be in the range between 
70 GeV and 105 GeV. 

Selection criteria on the energy of the most energetic electromagnetic cluster of the event and 
the acoplanarity and acollinearity of the lepton pair reject the residual Bhabha and radiative 
fermion pair-production backgrounds. 

The data and Monte Carlo distributions for the sum of Mu and M rec , the most discrimi- 
nating variable between signal and background, are displayed in Figure Eli. The yield of the 
selection for the total sample, and the separate energies are respectively given in Tables 2 and 3. 
The background is mainly constituted by non-resonant neutral-current four-fermion events. 

4 Results 

4.1 Measurement of the ZZ Cross Section 

The distributions of the variables presented in Figures [TJd, [Hi, Eb,Efc and Eli are separated into 
each centre-of-mass energy and are then fit to determine the cross section of the individual 
channels. 

A probability density function is built from the observed number of events in each of the 
bins of the distribution as a function of the signal cross section, fixing the background expec- 
tations. A flat positive distribution for its value is assumed. If a zero value of the cross section 
is contained in a ±34% confidence interval around the maximum of the probability density 
function, then a 95% confidence level upper limit is calculated. This maximum is otherwise 
quoted as the measurement, adopting this interval as the corresponding statistical uncertainty. 

Table 4 lists the results of all these fits together with the Standard Model predictions. 
Assuming these predictions as the relative weights of the different channels, the ZZ cross section 
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for each centre-of-mass energy can be calculated from a simultaneous fit to the five channels. 
The results of this fit are also presented in Table 4. All the measured cross sections agree with 
the Standard Model predictions. In the calculation of the cross section, the effect of the cross 
talk between the separate channels is found to be negligible. 
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0.01 
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0.79 
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0.04 


< 0.11 


0.01 


1.17 ±0.24 


0.92 


200 


< 0.24 


0.04 


0.06 ±0.04 


0.01 


1.25 ±0.25 


1.00 


202 


< 0.26 


0.04 


< 0.13 


0.01 


0.93 ±0.38 
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Table 4: Results, a* 1 *, of the individual and global cross section fits for the different centre-of- 
mass energies. The corresponding theory predictions, a th , are also given. Limits are at 95% 
confidence level. 

Figure presents the distribution of the reconstructed mass, M, for all the selected events, 
including those collected at lower centre-of-mass energies [HIIZI- A cut on the qqisi? and qqq'q' 
neural network outputs at 0.8 and 0.85 is applied, respectively. For the qq£ + £~ and the qqvv 
channels, M corresponds to M^c and Mfu, respectively. The average of the two di-jet masses 
is considered for the qqq'q' channel while for both the l + l~vv and £ + £~£' + l'~ channels the 
average of M« and M rec is used. A fit to the distribution of M is performed in terms of the 
ratio Rzz between the observed events and the predictions from Z pair-production and yields: 

R zz = 0.94 ±0.14 ±0.06, 

in agreement with the Standard Model. The first uncertainty is statistical and the second 
systematic, discussed in References |B] and and below. The cosine of the observed production 
polar angle 9 is presented in Figure for the same selected events. 

4.2 Study of Systematic Uncertainties 

Systematic uncertainties on the e + e~ — >ZZ cross section can be divided into sources correlated 
and uncorrelated among the channels. Their effects are estimated using the full 1999 data 
sample, and then propagated to the measurements performed at the different centre-of-mass 
energies. 

The main sources of correlated systematic uncertainty are the background cross sections and 
the energy scale of the detector. As they modify the shapes of the fit distributions, their effect 
is evaluated performing a new fit to calculate the cross section once their values are modified 
as listed in Table 5. Possible non-linearity effects for the energy scale are investigated. The 
effect of the uncertainty of the LEP beam energy is negligible. 

An uncertainty of 2% is attributed to the measured cross section to take into account the 
difference of the assumed relative weights of the different channels, given by the EXCALIBUR 
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calculation, with respect to those obtained with GRC4F, and to parametrise other uncertainties 
related to their calculation. 

Some sources of systematic uncertainty are uncorrelated among the channels but modify 
the shape of the output of the final neural network of the qqq'q' selection. These are the jet 
resolution, the charged track multiplicity and the b-tag, and their effect is presented in Table 5. 
The jet resolution includes a variation of ±2° on the jet direction and a smearing of t/34. A 
variation of the b-tag discriminant of ±2% models possible systematic effects and includes 
uncertainties in the Monte Carlo description of b-hadron jets. 



Systematic Source 


Variation 


5a zz (%) 


<fozz-bbx (%) 


Correlated sources 








WW cross section 


2% 


2.4 


2.6 


Four-jet rate 


5% 


2.1 
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Weu cross section 
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< 0.1 
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5% 
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Lep energy 
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see text 


0.3 
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1% 
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B-tag (qqq'q') 


2% 


1.6 


7.5 
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see text 


3.9 


3.3 


Simulation/Lepton Id 


see text 


2.5 


1.4 


Total 




7.5 


10.8 



Table 5: Systematic uncertainties on azz an d Czz->bbx- The total systematic uncertainty is the 
sum in quadrature of the different contributions. 
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0.7% 
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24.8% 
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1.5% 


1.3% 


0.7% 






Background MC statistics (^zz^bBx) 


8.1% 


0.6% 


1.3% 






Simulation/Lepton Id 


1.3% 


1.9% 


1.2% 


4.7% 


11.3% 



Table 6: Sources of uncorrelated systematic uncertainties on ozz and cr zz ^ bb x- 



Three additional sources of systematic uncertainty, uncorrelated among the channels, are 
considered: the Monte Carlo statistics of the signal and the background and the agreement 
between data and Monte Carlo. The latter comprises normalisation differences as derived 
from the comparison of data and Monte Carlo samples five to twenty times larger than the final 
ones, obtained by relaxing some selection criteria. It also includes differences in the shape of the 
distribution of the lepton identification variables around the adopted selection requirements. All 
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these uncertainties, listed in Table 6, do not affect the shape of the discriminating distributions 
and their effect on the total cross section propagates as summarised in Table 5. The total 
systematic uncertainty is the sum in quadrature of all these contributions. The measured cross 
sections then read: 

a zz (192 GeV) = 0.29 ± 0.22 (stat.) ±0.02 (syst.) pb (SM : 0.79 ± 0.02 pb), 
a zz (196 GeV) = 1.17 ± 0.24 (stat.) ±0.09 (syst.) pb (SM : 0.92 ± 0.02 pb), 
a zz (200 GeV) = 1.25 ± 0.25 (stat.) ±0.09 (syst.) pb (SM : 1.00 ± 0.02 pb), 
a zz (202 GeV) = 0.93 ±0.38 (stat.) ±0.07 (syst.) pb (SM : 1.03 ± 0.02 pb). 

The values in parentheses recall the Standard Model expectations. 

A new fit is performed in terms of the NC02 framework and the corresponding cross sections 
are derived as: 

cx£ z C02 (192 GeV) = 0.29 ±0.22 (stat.) ±0.02 (syst.) pb (SM : 0.77 ± 0.02 pb), 
a^ 002 (196 GeV) = 1.18 ± 0.24 (stat.) ±0.09 (syst.) pb (SM : 0.90 ± 0.02 pb), 
a£ z C02 (200 GeV) = 1.25 ± 0.25 (stat.) ±0.09 (syst.) pb (SM : 0.98 ± 0.02 pb), 
<t£ z C02 (202 GeV) = 0.95 ± 0.38 (stat.) ±0.07 (syst.) pb (SM : 1.01 ± 0.02 pb). 

The Standard Model expectations given in parentheses are calculated with the ZZTO [2*o] 
program and are assigned a ±2% uncertainty j2E|- The YFSZZ [22j package yields compatible 
estimations. As the relative weights of the different final states are set according to the Z boson 
branching fractions into fermions, the systematic uncertainty no longer includes the ±2% due 
to their predictions. On the other hand, a ±2% uncertainty is assigned to account for possible 
effects due to the extrapolation to the NC02 framework of the efficiencies and background 
estimations from the Monte Carlo simulations described above. 



4.3 b Quark Content in ZZ Events 

Z pair-production with at least a Z decaying into a b quark pair constitutes an interesting test 
bench of the detector capabilities to observe the minimal or a supersymmetric Higgs boson. 
These would in fact be seen as events with heavy particles decaying into a b quark pair, recoiling 
against a Z boson. Moreover, for the Higgs mass ranges under current investigation at LEP, 
the cross sections of these processes are similar. 

The qqq'q' final state analysis presents a high sensitivity to final states containing b quarks, 
as shown in Figure |2t> for the bbqq response of the neural network used to select the qqq'q' 
final states. 

The qqvv and qq£ + £~ selections, as summarised by the distributions of M 5 c for qq£ + £~ 
and the neural network output for qqvv, are complemented with the same b-tag variable as 
the qqq'q' selection. Its value is recorded for selected data and Monte Carlo events for each of 
the two hadronic jets. A single discriminant is then built for each channel from its two b-tag 
variables and the selection one. First the variables are mapped to achieve uniform distributions 
for the background. Then the product of their observed values is calculated event by event. Fi- 
nally the confidence level is calculated for the product of three uniformly distributed quantities 
to be less than the observed product. This confidence level is expected to be low for signal and 
flat for background. The final discriminant is the negative logarithm of this confidence level 
and is presented in Figure |U 
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The cross section calculation for the individual channels is performed as for the inclusive 
modes, considering the distributions in Figures Eb and|H and yields the results listed in Table 7. 
The combined result for 0"zz-«-bbx i s: 

^zz^bBx(192 - 202 GeV) = 0.31 ± 0.09 (stat.) ± 0.03 (syst.) pb, 

in agreement with the Standard Model expectation of 0.27±0.01 pb. In all fits, the contribution 
from other Z pair final states are fixed to their Standard Model expectations. The systematic 
uncertainties are evaluated in the same way as for the total cross section, and are presented in 
Tables 5 and 6. 







bbuu 


qqbb 


Measured cross section (pb) 
Expected cross section (pb) 


0.07 ±0.04 
0.031 


< 0.08 
0.057 


0.26 ±0.07 
0.178 



Table 7: Results of the individual ZZ- 
level. 



bbX cross section fits. The limit is at 95% confidence 



Figure El displays the measured total and bbX cross sections and their expected evolution 
with y/s, including data at lower centre-of-mass energies [HJ[7] and the theory uncertainties 
discussed above. 
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Figure 1: Distributions for data and Monte Carlo of: a) Invariant mass Mm of the lepton pair 
for the qq£ + £~ final state before the application of the cuts indicated by the arrows, b) The 
fit mass M 5 q of the qq£ + £~ final state, c) The mass Mf it of the hadronic system of the qquu 
final state after a kinematic fit that imposes the Z mass to the event missing four-momentum, 
d) The output NNo u t of the qqvv neural network. 
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Figure 2: Distributions for the qqq'q' selection of the outputs a) NNo u t4J °f the first neural 
network; the ZZ signal is superimposed with a cross section five times larger than the predicted 
one and the arrow shows the cut, b) NNoutZZ of the final neural network; signal expectations 
for events with no or at least one b quark pair are presented separately. Distributions of the 
sum of the visible and recoil masses for c) the £ + £~i/u and d) the £ + £~£' + £'~ selections. Data 
and Monte Carlo are shown. 1 a 
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Figure 3: Distributions in data and Monte Carlo at all the LEP centre-of-mass energies above 
the Z pair production threshold of a) the reconstructed mass M and b) the cosine of the 
production angle 9. Cuts on the qqw and qqq'q' neural network outputs are applied as 0.8 
and 0.85, respectively. 

17 



cd" qq 

o £ 
g- CD 

h3 



o 



o 

l-S 



p 3 

CO p. 

y. a 



p" 
cr 

o 

CD p 

| & 
!zb P 

9 P 



o 
a 



l-S 

o 



O 
p 

i-i 

O 



CD 

cr 
— i 

+ 

I 

P 

a 



L3 Vs=1 92-202 GeV a) 



i 



10 



LO 

■ 

o 

H— > 

c 

> 
LU 



Data 

ZZ^bbll MC 
Back. MC 



2 4 
Discriminant 



6 



tr 

o 

— 

en 



L3 Vs=1 92-202 GeV b) 



Data 

ZZ^bbvv MC 
Back. MC 



4i 



^1 



2 4 6 

Discriminant 



s 1 - 5 

c 
o 

O -I 

a> 1 

tf) 

</> 
</> 

o 



o 



0.5 



ZZ Data 

e e ->ZZ 
ZZ^bbX Data 
ee "^ZZ^bbX 




170 



180 190 

Vs (GeV) 



200 



210 



Figure 5: Measurements of the e + e — >ZZ and ZZ — > bbX cross sections, where statistical and 
systematic uncertainties are combined in quadrature. 
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